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Preferential attachment is a central paradigm in the theory of complex networks. In this contribution we consider 
various generalizations of preferential attachment including for example node removal and edge rewiring. We 
demonstrate that generalized preferential attachment networks can undergo a topological phase transition. This 
transition separates networks having a power-law tail degree distribution from those with an exponential tail. 
The appearance of the phase transition is closely related to the breakdown of the continuous variable description 
of the network dynamics. 
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Complex networks are found in nature, socio-economic sys- 
tems, technical infrastructures, and countless other fields. The 
macroscopic properties of such networks emerge from the mi- 
croscopic interaction of many individual constituents. Vari- 
ous models of complex networks have been proposed 01-0] 
and statistical mechanics of complex networks has become an 
established branch of statistical physics. Since complex net- 
works are non-equilibrium systems they do not have to obey 
detailed balance and may show many fascinating features not 
found in equilibrium systems. 

In order to explain the power-law degree distribution of 
many complex networks Barabasi and Albert f3\ introduced a 
model that combines network growth and preferential attach- 
ment. Preferential attachment means that new nodes enter a 
network by attaching preferably to nodes having a high de- 
gree. The probability that a new node u establishes an edge to 
an old node v with degree fc,, is determined by a node attrac- 
tiveness function A^.; it is given by padd{v) =^/t„/Lw^A:„,- In 
the original Barabasi-Albert model 131 the attractiveness of a 
node is strictly-proportional to its degree, viz. ~ k. 

The linear preferential attachment function of the Barabasi- 
Albert model was introduced as an ad hoc ansatz without any 
fundamental justification and many generalizations are con- 
ceivable, including attractiveness functions with small pertur- 
bations. 



Ak=k + E{k,k*) with e(;t,r) e ^(1) , 



(1) 



random node removal, rewiring, etc. Here we will consider 
generalized preferential attachment models that include tun- 
able parameters that allow to adjust the relative importance 
of node removal, rewiring, etc. These parameters define the 
phase space of a model. 

We will show that generalizations of the Barabasi-Albert 
model as considered here can dramatically affect the degree 
distribution; it falls into one of two possible universality 
classes. At each point of the phase space the universality class 
is uniquely determined by a parameter j3 . For j5 >Q the tail of 
the degree distribution is given by a power-law, while if j3 < 
it is an exponential. The two universality classes are separated 
by a topological phase transition. The behavior at the phase 
transition ()3 = 0) depends on the specific model. 



We claim that our results are valid for a large class of gener- 
alized preferential attachment models, but let us consider two 
specific examples for illustrative purposes. 

Model I At each time step with probability r a randomly cho- 
sen node is removed from the network and with probabil- 
ity 1 a new node is added to the network by establishing 
m links to m old nodes. These nodes are chosen randomly 
by preferential attachment with the attractiveness function 
Q. The perturbation function may be any bounded func- 
tion e{k,k*) e such that Aj- > 0, e.g. 



e{k,k*) 



\k*k/m if k<m 
k* else. 



(2) 



For k>m function (|3 implements a constant shift of size 
k*. The network starts to grow at time f = with A^o nodes. 
We are interested in behaviour where the number of nodes 
is much larger than unity. For r = 1 the number of nodes is 
conserved so A^o ^ 1 . but for r < 1 the average number of 
nodes grows with time and one may start with A^o = ^(1)- 

Model II We start from a small initially given network. At 
each time step one of three different operations is per- 
formed with probability p, q, or r, respectively, such that 
p + q + r ^ I andO < p,q,r < l.||lj 

Edge insertion (probability p) There are m new edges in- 
serted. For each edge one end u is chosen randomly from 
the set of nodes without any preference, while the other 
end V is selected preferentially with probability propor- 
tional to A|(.|, = ky + 1 . 

Rewiring (probability q) m links are rewired by selecting 
a node u (with non-zero degree) and one of its neighbors 
V. The edge between u and v is removed and a new edge 
from M to v' is added, where v' is chosen with probability 
proportional to A,;. , = kyi + 1 . 

Node addition (probability r) A new node is added to the 
network by establishing m links to m old nodes v which 
are chosen with probability proportional to A<.,, = k^ + l. 

Let pk denote the degree distribution of a generalized pref- 
erential attachment model, i. e., the relative number of nodes 
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Figure 1 : Degree distribution for model I with ttie perturbation function 
1^1, node deletion rate r = I, and m = 2; sub-figure (a) exponential 
case, sub-figure (b) power-law case. 



having degree k averaged over many network instances. This 
distribution can be calculated by solving a master equation 
with the boundary values /:>_ i =0 and Pky„,,^+\ = 0, where A:max 
denotes the maximal possible degree in the network. The ther- 
modynamic limit corresponds to ^max — * °°- For model I the 
master equation reads l6|| 

Sk.m = -m^^Pk-\ + {ni^+rk+\^ Pk- r{k+l)pk+i , 

which contains the mean attractiveness 
A = Y.AkPk . 

k 



(3) 



(4) 



The simultaneous numerical solution of equations Q and 
gives either a degree distribution with a power-law tail or an 
exponential tail, as shown in Figure ^ For model II we can 
write down a similar master equation and its solution obeys 
either a power-law or an exponential tail, too, depending on 
the parameters p, q, r, and m. 



Analysis of generalized preferential network models in 
terms of a continuous variable approximation [TJ sheds light 
on the appearance of the phase transition. We will outline the 
generic mechanism that leads to the topological phase transi- 
tion first and we will consider specific models later on. 

In a continuous variable approximation degrees are treated 
as continuous variables and random quantities are represented 
by their means. The degree distribution /:><. is given by 0] 

-1 



Pk 



Nit) 



dk{i,t) 



(5a) 



Here N{t) denotes the size of the network at time t. Node 
addition and removal lead to an effective network model-de- 
pendent growth rate a, i.e., N{t) = at+No. The quantity 
D{i,t) takes into account random node removal at rate p. It is 
the probability that a node that entered the network at time i is 
still present at time t > i and is given by the initial condition 
D{i,i) = 1 and the equation 



dD{i,t) 
dt 



N{t) 



(6a) 



The function k{i,t) denotes the mean degree of a node at time 
t > i that entered the network at time / and has not yet disap- 
peared from the network. The time < f is given by the solu- 
tion of the equation k{ik,t) = k. The quantity \dk{i,t)/ dt\yjj^ 
can be interpreted as the number of nodes with degree k at 
time t (neglecting node deletion). 

For generalized preferential attachment networks the evo- 
lution of k{i,t) is given by an equation which has the generic 
form 



dkjij) _^km) , f[k{ij)] 
dt 



, •. ■ / N • (7a) 
Nit) Nit) 

The term likii,t)/Nit) is the hallmark of preferential attach- 
ment H, while f[kii,t)]/Nit) with /(x) e ^(1) appears as 
a consequence of perturbations of the strictly hnear attrac- 
tiveness function, edge rewiring, node deletion, or other pro- 
cesses. 

Before analyzing the equations (Bat . ( l6al . and dTal i in detail 
it is useful to introduce a new time scale t'. For a > we set 
t ~ exp(af') and in the case a = we set f = t'No. With this 
new time scale the number of node additions per unit-time is 
proportional to the network size. Now equations ( Bat . ( l6al l. 
and (l7al read (omitting primes) 



Pk 



Diik,t) 



exp[a(f 
dDii,t) 



-ik)] 



dkii,t) 



i='k 



dt 



-pD(/,f), 



dkii,t) 
dt 



= l5kii,t)- 
I 



^f[kjij)] 
II 



©) 
&) 



The explicit dependence on the network size has dropped out 
of equations i5h\ . i6h\ . and (17b > and we need not to distin- 
guish between growing (a > 0) and non-growing (a = 0) net- 
works. 
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The occurrence of the phase transition can be pinned down 
to a change in the solution structure of ilbl . We have to dis- 
tinguish the three cases j3 > 0, j3 = 0, and j8 < 0. First let 
us consider j3 > 0. Because f[k{i,t)] ^ (l) equation ( 17b > is 
dominated by term I. For k{i,t) ^ 1 the leading term of the 
solution of ( 17b > is given by 



k{i,t) ~exp[j3(f 



and there follows 



ln;t 



Ik 



Plugging (|8}, (|9}, and the solution of (I6b> . 

D{i,t)=exp[-p{t-i)], 



(8) 



(9) 



(10) 



into dSbt we find that the combination of exponentials with 
different exponents yields a power-law degree distribution 
pii ^ k^"^ with exponent 



p + a 



1. 



(11) 



The power-law exponent (II 1> is determined entirely by three 
rates, the network growth rate a, the node degree growth rate 
j3, and the node deletion rate p. It diverges as j3 ^ 0. 

It turns out that the exponential growth of a node's degree 
(|8} is the key ingredient for the emergence of a power-law 
degree distribution. For j3 = this ingredient is lost and the 
phase transition appears. The equation /3 = defines the crit- 
ical phase line that separates the two universality classes. At 
the phase transition term I in equation Mhl vanishes and as a 
consequence the qualitative behavior of the solution of ( 17b > is 
determined by term II and not universal. Because /(x) G ^ ( 1 ) 
function must grow sub-exponentially. The solution of 
Mb\ is given by the implicit equation 



'^'■.0 1 



Ak' : 



(12) 



For j3 < term I dominates equation ( I7bt again. But now 
the continuous variable approximation breaks down. In fact, 
the time (|9|l becomes greater than f; that means it is in the 
future. The continuous variable approximation is not able to 
predict the correct degree distribution for /3 < 0. 

Let us see how our generic continuous variable analysis 
translates to specific models. In model I in each time step 
a node is added to the network and with probability p = r a 
node is removed from the network, as a consequence its size 
grows at rate a = (1 — r). The function k{i,t) is given by 



dk{i,t) k{i,t) + e[k{i,t),k*] k{i,t) 



dt 



AN{t) 



N{t) 



(13) 



The first term takes into account node addition via pref- 
erential attachment with attractiveness function ([0 and the 
second term corresponds to random node deletion. After 
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Figure 2: Phase diagram of the degree distribution of model I with the 
perturbation function |^ and m — 3. Symbols indicate locations in the 
phase space for data shown in Figure|3l 



reordering equation (I13> and comparing with ( fTal we find 
/3 = [m/ (M) — l]r. For )3 > the degree distribution is given 
by a power-law with exponent y~\ /{[m/ [rA) — l]r} + 1. In 
order to determine the universality class of the degree distri- 
bution we have to calculate j3 and A by solving equations (|3j 
and (|4} numerically. If we choose the perturbation function (|2j 
then the phase space is spaned by r and k*, where the critical 
phase line is given by [m/{rA) — l]r = 0, see Figure |2] Di- 
rect simulations of the network dynamics show that in the uni- 
versality class of power-law degree distributions k{i^t) grows 
indeed exponentially, as predicted by (|8}, while in the expo- 



10" 



10 




6 
/ - / 



10 



12 



Figure 3: Mean degree of a node at time t > i (at rescaled time scale) 
for model I networks with perturbation function in = 3, and re- 
moval rates r = 0.5 (open symbols) and r = 1 (solid symbols), re- 
spectively. The diagram shows networks at the phase transition, in 
the exponential regime, as well as in the power-law regime; see Fig- 
urel^for the corresponding position in the phase diagram. Symbols 
show binned data obtained by a direct simulation of the network dy- 
namics, lines are given by analytic computations of k{i, t) at the phase 
transition fT^ . 
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Figure 4: Phase diagram of model II for m = 2. 

nential class k(i^t) grows sub-exponentially, see Figure|3] 

Moreover, our findings for model I put some known results 
for specific variations of tiie Barabasi-Albert model into a uni- 
fying framework. For a vanishing node deletion rate r = and 
e{k,k*) = k* the mean attractiveness equals A = 2m + k* and 
for any k* the deg ree distribution has a power-law tail with 
7= 3 + k*/m (lyi. In the special case e{k,k*) = with node 
deletion rate < r < 1 the mean attractiveness A equals the 
mean degree of the network, which is given by 2m/ (1 + r). It 
follows that the degree distribution is in the power-law class 
with7=(3-r)/(l-r)L6,ll. 

Model II can be analyzed in a similar way. In this model the 
network size grows at rate a^r=l— p — q, and the evolution 
equation for k{i,t) is given by 

dk{i,t) mp t7ip[k(i,t) + 1] 
dt ~N(t)^ N(t)Ck+l) 

mqk{i,t) mq[k{i,t) + l] mr[k{i,t) + l] 

N{t)k Nit){k+l) N{t)ik+l) ' ' 

where k = 2m{l — q)/r denotes the mean degree of the net- 
work. After some rearranging of equation \\A\ and comparing 
with ( I7al we identify 

1 (1 -p-q)[-{2m+\)q^+{Am+l~p)q-2m] 
^ ~ 2 (l-^)[(2m+l)^-(2m + l)+p] ' 

So, for j3 > the degree distribution will follow a power-law 
with exponent 

{l-ci)[{2m+l)q-{2m+\)+p] ^ ^ 
— {2m + 1 + [Am + I — p)q — 2m 

and the critical phase line is determined by the equation 
(2m + 1 )q^ — (4ot +1 — p)q + 2m = 0, as shown in Figure^] 



To summarize, we have investigated a broad class of gen- 
eralized preferential attachment models and shown that small 
perturbations of the attractiveness function Q, node removal, 
or edge rewiring cannot only change thepower-law exponent 
of the original Barabasi-Albert model |5|]. These variations 
can even lead to a qualitative change of the degree distribu- 
tion from a power-law to an exponential. 

It had been argued 1^ "that net growth is one of the funda- 
mental requirements for the generation of power-law degree 
distributions" of preferential attachment models. However, 
we have shown, that it is not the growth of the network it- 
self that leads to a power-law degree distribution. The mean 
degree of individual nodes has to grow sufficiently strongly, 
i. e. j3 > 0, to generate a power-law degree distribution. 

We also demonstrated that negative perturbations of the lin- 
ear attractiveness function as in model I have a tendency to 
stabilize power-laws. In the light of this fact and the ubiquity 
of power-laws in nature it would be interesting to know if in 
real world networks such negative perturbations do exist. 
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